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ASTEROSEISMIC SIGNATURES OF SMALL CONVECTIVE CORES 
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ABSTRACT 

We present an analytical study of the effect of small convective cores on the oscillations of solar-Uke pulsators. 
Based on an asymptotic analysis of the wave equation near the center of the star, we derive an expression 
for the perturbations to the frequencies of radial modes generated by a convective core and discuss how these 
perturbations depend on the properties of the core. Moreover, we propose a diagnostic tool to isolate the 
predicted signature of the core, constructed from a particular combination of the oscillation frequencies, and we 
validate this tool with simulated data. We also show that the proposed tool can be applied to the pulsation data 
soon expected from satellite missions such as CoRoT and Kepler to constrain the amplitude of the discontinuity 
in the sound speed at the edge of the convective core, the ratio between the sound speed and the radius at this 
same location, and the stellar age. 

Subject headings: convection — methods: analytical — stars: interiors — stars: oscillations 



1. INTRODUCTION 

Space-based asteroseismology began in 1999, when the 

5 cm star tracker mounted on the Wide-field InfraRed Ex- 
plorer (WIRE) s atellite was used to detect oscillations in the 
K giant a UMa (iBuzasi et al.ll200dl) . and this instrument has 
continued to produce a stea dy stream of asteroseismic data 
(e.g., see iBruntt et al.l l2007h . The 15 cm Micro-variability 

6 Oscillations of Stars (MOST) satellite was launched in 
2003, and was the first satellite in orbit that was actually de- 
signed for asteroseismology (Wal ker et al.ll2003h . Although 
MOST's failure to de tect solar-like oscillations in Procyon 
(iMatthews et al.l l2004l) was later traced to a lar ger than ex- 
pected source of non-white instrumental noise (B edding et al.l 
|2()05), its ability to conduct nearly uninterrupted time-series 
photometry for durations of up to two months has produced 
exquisite data for pulsating stars acro ss the H-R diagram (e.g. , 
see Randall et al. 2005; Aerts et all 120061; ISaio et al.1 120061; 
iRowe et al.ll2006t ICameron et al.ii2007l) ! 

The next era of space-based asteroseismology began 
with the successful launch of the 27 cm Convection, Ro- 
tation and plan etary Transits (CoRoT) satellite in Decem- 
ber 2006 (Bagl in etaLll2006h . and will continue with the 
anticipated launch of the 95 cm Kepler satellite in 2008 
(IChristensen-Dalsgaard et al. 2007). These missions promise 
advances in the study of solar-like oscillations comparable to 
what WIRE and MOST have achieved for the larger amplitude 
classical pulsators. With the frequency precision expected to 
approach ^0. 1 /iHz from this next generation of satellites, we 
can begin to consider the possible detection of extremely sub- 
tle asteroseismic signals. 

There is a long history of constructing combinations of the 
observed oscillation frequencies to extract useful asteroseis- 
mic information. The simplest such combination is the aver- 
age frequency spacing between modes of consecutive radial 
order but the same spherical degree — the large separation — 
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which provides an estimate of the mean stellar density dCoxl 
IT980I) . Similarly, the average spacing between modes of con- 
secutive radial order but with spherical degrees that differ 
by 2 — the small separation — is sensitive to chemical gradi- 
ents in the d eep interior, tracing the stellar evolutionary state 
(iDappen et al. 1988). 

The idea that sharp variations in the internal structure of 
distant stars generate signatures in the frequencies of low de- 
gree modes that can be recognized and interpreted with the 
appropriate seismic tools has been explored in a number of pa- 
pers. In this context, theoretical derivations were carried out 
of the signal expected from sharp v ariations taking place at the 
boundary of convective envelopes ("Monteiro et al."2000) and 
in the region of hehum ioniz ation (Monteiro & Thompso^ 
119981; iHoudek & Goughll200'7l) . Moreover, the expected sig- 
nature of convective cores on the small frequency sep arations 
has also been considered by [Roxburgh & Vorontsovl (l200ll) . 
However, all of these works assumed that the sharp variation 
in the structure under investigation was located far from the 
turning points of the oscillation modes, i.e., well within their 
propagation cavities. Consequently, the analyses carried out 
in these works do not directly apply to studies of small con- 
vective cores, such as those present in main-sequence solar- 
like pulsators. 

Attempts to identify and interpret in a systematic manner 
the signatures of convective cores in solar-like pulsators di- 
rectly from simu l ated data have recently been carried out by 
iMazumdar et al.l ( l2006i) . Based on an analysis of the results 
obtained from simulated data, the authors suggested a diag- 
nostic tool to estimate the masses of the convective cores 
and the stellar ages. While recognizing the importance of 
studying simulated data, the chances of constructing the op- 
timal diagnostic tool to detect and characterize small con- 
vective cores may be greater if we know what signal to ex- 
pect. Moreover, a theoretical understanding of the expected 
signal can provide us with a direct link between the struc- 
ture under investigation and the observations — moving us 
toward our goal of finding the best way to infer informa- 
tion about that structure. With this in mind, in this paper 
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we derive the expected signal of a small convective core on 
the oscillation freque n cies of solar-hke pulsators. U nlike 
Monteiro et al.' (200(f), ' Roxburgh & Vorontsovl (1200 lb . and 
Houdek & Gough (2007), we do not assume that the region 
of sharp structural variation, which in the present case is the 
edge of the convective core, is placed well within the propa- 
gation region of the modes. In fact, the analysis is focused on 
the opposite case, in which the edge of the convective core is 
sufficiently close to the center to affect the radial modes more 
significantly than the modes with degree £>1. 

We derive the theoretical expression for the signature of a 
small convective core in ^2.11 In 32.2l we construct a diagnos- 
tic tool, based on a particular combination of the frequencies, 
to isolate this signature. In ^we validate the theoretical pre- 
dictions and the diagnostic tool using simulated data from a 
series of 1.3 M0 models. Moreover, we discuss the successes 
and limitations of the proposed tool, as well as the informa- 
tion that it can reveal. Finally, in ^we outline the potential 
for future work in this area, and we discuss possible observa- 
tional tests. 

2. SIGNATURE OF A CONVECTIVE CORE 

2.1. Theoretical expectation 

Let ^ = (^,0,0)re""^' be the displacement, defined with re- 
spect to a spherical coordinate system (r, 6*, 0), for radial, adi- 
abatic, acoustic oscillations in a spherically symmetric star, 
with ijj the oscillation frequency and t the time. Through 
the appropriate mathematical manipulation of the equations 
of motion and mass conservation for radial pulsations, under 
the adiabatic approximation, the dimensionless displacement 
£ is found to obey the second order differential equation (e.g. 
lGoughllT993h . 



e = o, (1) 



where 7 is the first adiabatic exponent, p is the pressure and 
p is the density. By defining new dependent and independent 

variables, S = rH^r'l^ p^l^cS, and z = ln(r/ri-ef) respectively, 

where c = \fjplp is the sound speed and rret is a fiducial 
value of r, and substituting into Eq. ([T]i, the latter reduces to 
the standard form. 



d^S 
dz2 



+/:2s = o, 



where 



1 

■4' 



(2) 



(3) 



and (jjc is a critical acoustic frequency (e.g., see iGoughll 19931 
Eq. 4.8.7). Since our aim is to study the effect of a small 
convective core on the eigenfrequencies, we are particularly 
interested in the form of u>c close to the center of the star. 
From the general expression for lJc, we find that when r — > 0, 
K, Ic^ jp-. Using this expression for cj^, we then find the 
critical frequency w* at which = 0\a be. 



3c 
Yr' 



(4) 



A model with a small convective core exhibits a sharp struc- 
tural variation close to the inner turning point of Eq. This 
sharp variation can be seen in the sound speed profiles of 
main-sequence models with mass M = L3 Mq, shown in Fig- 
ure [T] for a range of stellar ages. Our aim is to determine 
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Fig. 1. — Sound speed in the inner layers of main-sequence models with 
M = 1 .3 Mq , and ages ranging from 0.25 Gyr (top curve) to 5.25 Gyr (bottom 
curve) in intervals of 0.5 Gyr The discontinuity in c marks the edge of the 
convective core, which expands during the evolution of the star up to an age 
of 3.25 Gyr and then begins to contract. 



how the eigenfrequencies of an otherwise similar model with 
a smooth structure around the edge of the convective core — 
hereafter the unperturbed model — are modified by the pres- 
ence of this sharp variation. 

To achieve this we must find an expression for S that is valid 
close to the inner turning point of Eq. Such an expression 
for S is derived in Appendix lAl In short, we define new de- 
pendent and independent variables and substitute them into 
Eq. ^ to obtain an equation that resembles the Airy equation 
with an additional term (cf. Eq. lA3l in AppendixlA]). Applying 
Olver's comparison method ( 01ver..l974.) we then find, close 
to the inner turning point. 



S«fl|x|'/4|7i:|-i/2[^.(0)_A'i(0)jc] 



(5) 



where a is a constant related to the normalization of eigen- 
functions, Ai(0) and A'i(O) are the values taken at x = of the 
first solution to the Airy equation and its first derivative re- 
spectively, and X is the new independent variable, which is 
related to z through Eq. (lAlb in Appendix [A] In deriving 
Eq. (|5]l we have neglected terms of (9(x^). So, this equation is 
approximately correct when x is sufficiently close to 0. Note 
that unlike r, the independent variable x depends on the value 
of the eigenfrequency. For a given frequency it goes through 
at the inner turning point, and is negative when K} <0 and 
positive when > 0. The relation between x and r for dif- 
ferent values of the eigenfrequency uj, as well as the regions 
where <0 and > 0, are shown in Figure |2] for an un- 
perturbed stellar model with M = 1.3 Mq. The unperturbed 
models were obtained by fitting the sound speed profiles of 
each stellar model with smooth functions around the edge of 
the convective core. The density profile was then derived by 
assuming hydrostatic equilibrium. 

To study the effect of a small convective core on the oscil- 
lation frequencies, we start from an integral equation for the 
frequency, which can be derived from Eq. ([T]) after multiply- 
ing by ^ and integrating by parts between r = and r = R, 
where R is the radius at the stellar surface. The resulting well 
known equation is, 

(6) 
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where 



h = 



dr' 



dr, 



and 



B = B{R)-B{Q)-- 



dr 



r=R 



r=0 



When the oscillations are below the acoustic cutoff frequency, 
the surface term B can safely be neglected. In fact, it follows 
from the regularity condition at r = that B{Q) = and, given 
the small value of p{R), if the oscillations are evanescent in 
the outer layers B{R) will also be small. 

We now consider the above mentioned unperturbed smooth 
model and an otherwise similar model with a sharp structural 
variation at the edge of a small convective core — hereafter the 
perturbed model. If B is neglected in Eq. (|6]l, the difference 
between the oscillation frequencies in the perturbed and un- 
perturbed models in a linear approximation is given by 



5lo ■■ 



6h 



(7) 



where all unperturbed quantities refer to the smooth model. 
Equation ^ with B neglected c onstitutes a va riational prin- 
ciple for the frequency uj (e.g. iChandrasekhaiil 1 9641: iGoughl 
Il993h . Thus, to first order the perturbation to the eigenfunc- 
tions does not contribute to the perturbation to the eigenfre- 
quencies and we can write. 



and 



dh= / Spr'^^^dr, 
Jo 



(8) 



r'-[3S(^p)-4Sp]e}dr, 



(9) 



where Sp, Sp and Sj denote perturbations at fixed radius. 
To obtain an expression for Siv as a function of w, we ma- 
nipulate Eqs. (|7]|9]l in a way very similar to that presented 
in Monteiro et al. (1994). In particular, Eqs. dHJ and (|9]l are 
written in terms of the perturbations (5(7/7), Sc^ and Sg, with 
g = —-^^ and then integrated by parts (see Appendix FBI for 
details). As discussed in Appendix iBl the perturbation to the 
eigenfrequencies is determined predominantly by the terms 
associated with Sc^. Describing thi s pert urbation with a mod- 
ified step function, we find (cf. Eq. lBlOb 
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(10) 



where As is a positive constant related to the sound speed in- 
crease at the edge of the convective core and is the value 
taken by the independent variable x at the location of the edge 
of the convective core. Both x and Xd are functions of fre- 
quency. Through this dependence, Eq. ( fTOl i expresses how the 
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Fig. 2. — The full curve shows the critical frequency obtained by im- 
posing K- = for the inner layers of an unperturbed main-sequence model 
with M = 1.3Mq at the age when its convective core reaches maximum 
size (3.25 Gyr). For each frequency, this curve gives the maximum propa- 
gation depth of a radial mode when the eigenfunctions are the solutions to 
Eq. (2). The vertical line shows the radial position of the edge of the con- 
vective core. The horizontal lines are axes for the independent variable x in 
the range [—0.5, 0.5]. These are shown for eigenmode frequencies of the cor- 
responding perturbed model, with radial orders ranging from « = 1 1 (bottom 
line) to n = 20 (top line). 



frequency perturbation is modulated with frequency. Accord- 
ing to this equation the frequency perturbation goes through 
zero when .^d = 0, i.e. when the frequency is such that the turn- 
ing point and the edge of the convective core are located at the 
same depth. 

If the signal anticipated by Eq. ( fTOl ) could be detected in the 
oscillation frequencies of real stars, then information about 
the size and position of the sound speed discontinuity at the 
edge of the convective core could be inferred from the ampli- 
tude As and the frequency at which 6uj = 0, respectively. With 
this in mind, in the next section we present a diagnostic tool 
that will allow us to isolate the signal derived above. 

2.2. Diagnostic tool 

When deriving Eq. ( fTOl i we assumed that our model dif- 
fers from a smooth model only at the edge of the convec- 
tive core. However, sharp variations in the internal structure 
of a star are known to occur also at other locations, such as 
the base of the outer convective region and the region of he- 
lium ionization. These, in turn, will generate additional per- 
turbations to the eigenfrequencies, as compared with those 
of a smooth model ( e.g., Monteiro et al. 1994; Basu 199*^ 
iMonteiro & Thompso n 2005: Houdek & Goughl2007h . Thus, 
to isolate the signal produced by the edge of a convective core, 
we must combine the oscillation frequencies of real stars in 
such a way as to cancel out, as much as possible, the signal 
coming from all other sources. With this in mind, we write 
the frequency v„ i = LOn^/ln of a mode with radial order n 
and degree ^ in a given star as the sum of three components 



(11) 



where v^^ ^ is the frequency of the mode in the smooth model, 
5v'^ ^ is the perturbation to that frequency arising from all 
sharp variations in the structure of the star taking place well 
within the propagation region of the same mode (which we 
shall consider to be of low degree, with I < 3), and Si^^ ^ is 
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the frequency perturbation produced by the edge of the small 
convective core. 

While sharp variations in the envelope of a star will affect 
the frequencies of all modes of low degree, the discontinu- 
ity in the sound speed at the edge of a small convective core 
for a star slightly more massive than the Sun might be ex- 
pected to affect primarily the frequencies of radial modes. 
This is not to say that the frequencies of modes of degree 
£= I, which also propagate to very deep layers, will not be 
affected at all by the convective core. In fact, to correctly 
derive the expected signal for a dipolar mode, the Cowling 
approximation should be avoided. The derivation shou ld then 
be started from the formalism recently presented by iTakatal 
(|2006), which shows that the inner turning point of dipolar 
modes is modified when the equations are treated without ne- 
glecting the perturbation to the gravitational potential. How- 
ever, such an analysis is beyond the scope of the present pa- 
per Hence, based on the idea that radial modes are most 
sensitive to the deepest layers of a star, we will proceed by 
assuming Suf^^ ~ for modes with £ > I. This assumption 
will be tested in the following section. Moreover, given that 
we are interested in modes with degree £ < 3, we neglect the 
weak dependence of Si^j'^^ on £ (e.g.. i Monteiro et al] Il994t 
iHoudek & G ough 2007) when subtracting the frequencies of 
nearly degenerate modes. With these assumptions, the scaled 
small separations = — i^n-i^+a) / (4£+6) for pairs 

of modes with degree € < 3 divided by the large separations 
Az^„.£ = Vn+\.i-Vni take the form 



Do 



v. 



«,o 



- v. 



«-l,2 



5k 



«,o 



£>13 



(12) 
(13) 







In the high-frequency asymptotic regime, the quantity 

^^<iZ<iM^ (14) 
4£+6 Ai/„,£ 

in the unperturbed model is essentially independent of 
the pair of modes considere d (£ = 0, 2 or £ = 1,3; see 
[Roxburgh & VorontsovllIOOObllal l2003h . Thus, we proceed by 
subtracting Eq. ( fT3T l from Eq. (fT2l) to find the relation 



Dq2 Di 



5v: 



n.O 



Aj/„_i,i Aj/,,,0 6Ai/„_i.i 



(15) 



Note that for the large separations appearing in the denomi- 
nators, we used pairs of modes that encompass the frequency 
range of the modes used in the scaled small separations, which 
appear in the numerators. When defined in this way, the 
frequency differences presented in Eq. ( fTsT i are a smoother 
function of frequency than when Ai^,, o and /S.v„ \ are used in 
the denominators of the first and second term, respectively. 
In practice, this means that when subtracting Eq. (flJl l from 
Eq. ( fT2] i we also had to assume that f^'^"" w 

The frequency perturbation Sv^ g is expected to be modu- 
lated with frequency according to Eq. ( fTOb . Thus, if the as- 
sumptions leading to Eq. (fTsT i hold, we would expect to see 
a corresponding modulation when the oscillation frequencies 
observed in a given solar-like pulsator are combined in the 
manner prescribed. In the next section we test this predic- 
tion on simulated data and we discuss what might be inferred 
about solar-like stars from this diagnostic tool, once observa- 
tions with the necessary precision become available. 
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Fig. 3. — Frequency differences, as defined by the left hand side of Eq. )15t . 
for a model with M = i.3 Mq at an age of 0.25 Gyr when its convective core 
is very small and there is no discontinuity in the sound speed (see the top 
curve in Figure [T). 



3. TESTS WITH SIMULATED DATA 
3.1. Models 

We used the Aarhus ste llar evolution code (ASTEC; 
IChristensen-Dalsgaardll 1 982l) coupled with the adiabatic pul- 
sation code (ADIPLS) to calculate the theoretical oscilla- 
tion frequencies for 1.3 Mq models with ages from 0.25 
to 5.25 Gyr, sampling the full range of convective core 
sizes along the evolution tra ck. The ASTE C code used the 
equation of state (EOS) of lEggleton et alj (Il973h without 
Coulomb corrections, and opacities from the OPAL tables 
(Iglesias & Rogers 1996), supplemented by Kurucz opacities 
at low temperatures. The nuc lear reaction rates came from 
iBahcall & Pinsonneaultl (Il992h. convection was d escribed by 
the mixing-length theory o f lBohm-Vitens^ (119581) . and we did 
not include the effects of diffusion. 

We set the initial metallicity to Zq = 0.02, and the initial 
hydrogen mass fraction to = 0.74. We fixed the mixing- 
length parameter at a = 1.9 and included core overshoot us- 
ing aov = 0.25 with complete mixing in the overshoot region 
and considering the associated changes to the adiabatic tem- 
perature gradient. For each selected model, we calculated 
the radial and non-radial p-mode frequencies for oscillations 
with spherical degree ^=0-3 and radial order n=1^30. We 
used these frequencies to construct the combinations shown 
in Eq. ( fTsT l. and compared them to the theoretical predictions. 

3.2. Testing the diagnostic tool 

A number of approximations were made when deriving the 
diagnostic tool proposed in Eq. ( fTSl l. Thus, it is essential that 
we evaluate the ability of this seismic tool when it comes to 
isolating the signal from the convective core. 

With this in mind, we start by testing the last assumption 
made when deriving Eq. ( fTSl ) — namely, the weak dependence 
of the value of Ds in the unperturbed model on the pair of 
modes considered. To this end, we compute eigenfrequencies 
for our youngest model. At an age of 0.25 Gyr the convective 
core is extremely small, and there is no apparent discontinuity 
in the sound speed (see the top curve in Figure[T]i. Under these 
conditions, we expect 5v'^^ o ~ 0- Consequently, if the quantity 
Ds were strictly independent of the pair of modes considered, 
the frequency differences defined by the expression on the left 
hand side (LHS) of Eq. ( fTSl l would be approximately zero for 
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Fig. 4. — Left panel: the same as Figure|3] but for models with ages ranging from 0.25 Gyr (bottom left) to 5.25 Gyr (top left) in intervals of 0.5 Gyr. En'or bars 
are shown for the model with maximum convective core size. These correspond to Icr errors in the frequency differences when the relative error in the individual 
frequencies is assumed to be 10~*. For this model, the radial order of the modes spanning the frequency range considered varies between n = 9 and n = 27. Right 
panel: comparison between the theoretical signal expected from the convective core (dashed lines) and that contained in the frequency differences (full curves) 
for models with ages ranging from 2.25 Gyr (bottom left) to 3.75 Gyr (top left) in intervals of 0.5 Gyr. To determine the expected signal, we computed 5v'^^ ^ 
using Eq. jlOt and then normalized it by the large separations to reproduce the right hand side of Eq. )15t . 



all frequencies. 

Figure |3] shows the frequency differences obtained for our 
youngest model. Clearly, in the lower frequency domain the 
frequency differences deviate significantly from zero. This 
reflects the fact that Dg is nearly independent of the pair 
of modes considered only in the high-frequency asymptotic 
regime. At lower frequencies this assumption breaks down, 
and consequently so does the approximation made when de- 
riving Eq. ( fTSl l from Eqs. (fT2l l and ( fT3] l. Nevertheless, at high 
frequencies (> 2 mHz in Figure |3]l the frequency differences 
converge to a constant small value, which we shall designate 
I'asymp, Confirming our assumption of the near degeneracy of 
Dg with respect to the pairs of modes ^ = 0,2 and t= 1,3. 

To compare these results with those obtained for older mod- 
els, we determine the frequency differences defined by the 
LHS of Eq. ( fTSl ) for our sequence of models, with ages rang- 
ing from 0.25 Gyr to 5.25 Gyr The results are shown in 
the left panel of Figure |4] It is clear from this figure that in 
the high-frequency domain, where Eq. (flST l is expected to be 
valid, the absolute value of the slope of the frequency differ- 
ences increases with age. This can be understood as a direct 
consequence of the increase with age of the sound speed dis- 
continuity at the edge of the convective core. According to 
Eq. ( [Tol l, such an increase manifests itself through the ampli- 
tude As by an increase in the absolute value of the frequency 
perturbation at fixed frequency. This increase is consequently 
reflected in the frequency differences defined by the LHS of 
Eq. ( fTsT i, and thus in the slope of the curves plotted in Fig- 
ure |4] 

To determine whether the signature of the convective core 
on the frequency differences might be detected in future ob- 
servations of solar-like pulsators, we have included for one of 
our models the expected Icr error bars for the frequency dif- 
ferences, calculated with the assumption that the relative error 
on the individual frequenci es is 10"^ (^■^■' 0-^ mHz at 2 mHz 
as expected from CoRoT; iBagUn et al.ll2006h . It is evident 
from the figure that if individual frequencies are indeed deter- 
mined with such precision, we will be able to not only detect 
the effect of the convective core on the frequencies, but also 
distinguish between different possible models based on this 



signature. Moreover, we see that our ability to detect this sig- 
nature in real data will depend on the number of radial modes 
detected, as well as on their radial order. 

We made additional assumptions, other than the one re- 
ferred to at the beginning of this section, while deriving 
Eq. ( fTSl l. In particular, we assumed that the effect of the con- 
vective core on the frequencies of modes with degree l>\ 
is negligible, when compared to the effect on the frequencies 
of radial modes. To check whether this assumption is valid, 
we compare the frequency differences calculated for a subset 
of our models with the theoretical expectations derived from 
Eq. ([Tol l. The results are shown in the right panel of Figure|4] 
We see from the figure that in the high-frequency asymptotic 
regime the theoretical signal derived from the analysis of the 
wave equation near the inner turning point is indeed isolated 
when the oscillation frequencies are combined in the way sug- 
gested by Eq. (flSl l. We note that there is no model dependent 
free parameter in the theoretical expression for the frequency 
perturbations. The difference in the slopes of the theoretical 
curves at high-frequency, seen also in the frequency differ- 
ences calculated from simulated data, is thus a direct conse- 
quence of the dependence of the amplitude A ^ on the increase 
in the sound speed at the edge of the convective core. 

In contrast to these results. Figure H] shows that in the low- 
frequency domain the combination of frequencies suggested 
by Eq. ( fTSl l does not properly isolate the theoretical signa- 
ture expected from Eq. (fTol i. Just as in the case discussed 
for Figure [3] this is likely a consequence of the breakdown 
at low frequencies of the assumption that Dg is independent 
of the pair of modes considered. According to Figure[3] when 
the frequency perturbation associated with the convective core 
is negligible, the frequency differences become negative in 
the low-frequency domain, decreasing as the frequency de- 
creases. On the other hand, the frequency perturbation asso- 
ciated with the convective core becomes positive for frequen- 
cies below a particular value — namely, the frequency at which 
the position of the discontinuity in the sound speed coincides 
with the inner turning point of Eq. (|2]i. In the low-frequency 
domain, we thus expect our diagnostic tool to show a superpo- 
sition of these two behaviors. Indeed, this is seen for frequen- 
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Fig. 5. — Left panel: measurement of the slope of the frequency differences in the high-frequency domain as a function of the relative increase in the squared 
sound speed at the edge of the convective core (see text for the mathematical definition of A). The ages of the models shown range from 2.25 Gyr (top left) to 
5.25 Gyr (bottom right) in intervals of 0.25 Gyr Right panel: for models with ages ranging from 2.25 Gyr to 5.25 Gyr, we show the frequency j^max at which the 
frequency differences suggested by the LHS of Eq. U5t reach their maximum (crosses), and for the same models the value of (27r)"' 3c/2r (circles) at the edge 
of the convective core. Since c is discontinuous at the edge of the convective core, we use the mean of the values of c on the left and right side of the discontinuity 
when calculating this expression. For models with ages greater than 4 Gyr, we show the frequency i/cross at which the frequency differences cross the asymptotic 
value for the youngest model (diamonds). 



cies smaller than ^-^1.5 mHz in the left panel of Figure]?] The 
competition between these two behaviors results in the pres- 
ence of a maximum in the frequency differences calculated 
for each model. 

3.3. Characterizing the convective core 

The validation of the diagnostic tool proposed in Eq. ( ITSl l 
in the high-frequency domain makes us confident that we un- 
derstand the origin of the signature observed in the simulated 
data, and that we can use this signature to extract information 
about the convective core. As mentioned in ^2.11 and ^3.21 
information about the size of the sound speed discontinuity 
at the edge of the convective core can be extracted from the 
amplitude of the signal, and in principle information about 
the ratio [c/r] edge at the edge of the convective core can be ex- 
tracted from the value of the frequency at which the frequency 
perturbation goes through zero. Moreover, given the age de- 
pendence of the size of the sound speed discontinuity at the 
edge of the core, the signature discussed above might also be 
used to infer the evolutionary status of pulsating stars within 
the region of the H-R diagram considered here. 

These ideas are illustrated in Figure ]5] In the left panel 
the quantity A, defined as the difference between the fre- 
quency differences ^ ~ of modes with radial or- 
ders n = k + 2 and n = k-2, with k such that ^ 1750 /iHz 
for each model, is shown as a function of the relative increase 
in the squared sound speed at the edge of the convective core. 
The quantity A provides a measure of the slope of the fre- 
quency difference curves, such as those plotted in the left 
panel of Figure ]4] in the high-frequency domain. It is clear 
from Figure ]5] that a linear relation exists between A and the 
squared sound speed increase at the edge of the convective 
core. We also considered alternative quantities indicative of 
the slope of the same curves taken at constant radial order 
rather than at constant frequency. In all cases the results in- 
dicate a linear relation with the squared sound speed increase 
at the edge of the convective core. Note that when studying 
data with noise, the use of the quantity A defined by only two 
data points is not advisable. Instead, it would be more conve- 



nient to extract the slope of the frequency difference curves in 
the high-frequency domain through a least-squares fit to the 
appropriate part of the curve under consideration. 

Since the proposed diagnostic tool fails to reproduce the 
theoretical signal in the low-frequency domain, the idea of us- 
ing this tool to determine the value of the frequency at which 
the frequency perturbation associated with the convective core 
goes through zero — and thus extracting information about the 
physical size of the core — becomes compromised. Neverthe- 
less, since the maxima that characterize the curves defined 
by the frequency differences are a consequence of the pos- 
itive frequency perturbations in the low-frequency domain, 
we might expect the position of these maxima to contain in- 
formation about the size of the convective core. With this 
in mind, we have determined for each model the frequency 
i^m-dx at which the frequency differences defined by the LHS 
of Eq. (ITSl l reach their maximum value. These are plotted as 
a function of age in the right panel of Figure ]5] and are com- 
pared with the frequencies [(27r)"' 3c/2r] ^^^^ at which the fre- 
quency perturbation associated with the convective core goes 
through zero. 

Clearly, the frequency at which the frequency differences 
reach their maxima, indicated by the crosses in Figure ]5] is 
strongly correlated with the age of the model. This is not 
surprising since the magnitudes of the frequency perturba- 
tions, which influence the position of the maxima, are also 
correlated with age. Nevertheless, a change in the slope is 
noticeable around an age of ~4 Gyr This is likely a conse- 
quence of the contracting core. As seen in Figure H] when 
the core begins to contract, the sound speed variation near the 
edge of the convective core begins to spread over a larger ra- 
dial extent. Our theoretical analysis, which assumes that the 
sound speed variation at the edge of the core is well repre- 
sented by a modified step function, becomes less adequate in 
this case. Nevertheless, an inspection of the frequency dif- 
ferences shown in the left panel of Figure ]4] indicates that 
in the older models where this effect becomes evident, the 
amplitude of the frequency differences in the low-frequency 
domain is indeed substantially larger than in all other mod- 
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els. This may resuh from a sudden increase in the amplitude 
of the frequency perturbation in the low-frequency domain, 
or from an improvement in the degeneracy of Ds calculated 
for pairs of modes with degrees £ = 0,2 and 1= 1,3. What- 
ever the case, the relative importance of the latter in deter- 
mining the functional form of the frequency differences ap- 
pears to be smaller for models with ages greater than ^4 Gyr 
Thus, one might expect that for these models the frequency 
[(27r)"' 3c/2rJ^j^^ will be better estimated by the frequency 

i^cross at which the frequency differences cross the asymptotic 
value for the youngest model in Figure [3] Indeed this seems 
to be the case, as shown by the diamonds in the right panel of 
Figure m 

4. DISCUSSION 

We have derived the expected signature of small convective 
cores on the oscillation frequencies of solar-like pulsators, we 
have suggested a diagnostic tool to isolate this signature, and 
finally we have discussed what information the signature may 
reveal about the convective cores themselves and about the 
evolutionary status of the corresponding stars. 

The results demonstrate that our ability to isolate the con- 
vective core signature with the suggested diagnostic tool is 
significantly better in the high-frequency asymptotic domain 
than at lower frequencies. Consequently, the information in- 
ferred from this diagnostic tool at high frequencies is expected 
to be more robust. In particular, this is the case for the relation 
found between the slope of the frequency difference curves 
and the relative increase in the squared sound speed at the 
edge of the convective core. Based on the theoretical signal, 
which in the high-frequency domain is well represented by 
the proposed diagnostic tool, we expect this relation to hold 
also in a more general context when models of different mass, 
chemical composition, or with different amounts of core over- 
shoot are considered. 

By contrast, information inferred from the behaviour of the 
diagnostic tool at lower frequencies, such as the estimation 
of the ratio [c/r]edge from the position of the maximum in 
the frequency differences, or from the value of the frequency 
at which these differences cross I'asymp, is likely to show a 
stronger model dependence. As emphasized earlier, at lower 
frequencies the functional form of the diagnostic tool is influ- 
enced by the deviation from degeneracy of the quantity Ds- 
Consequently, this tool does not capture, as well as it does at 
higher frequencies, the effect of the convective core on the os- 
cillations, even though it is still influenced by it. Therefore, 
it is important to extend the present study to different models, 
varying quantities such as mass, chemical composition, and 



core overshoot, to determine the extent to which the quantities 
studied — ^particularly in the low-frequency domain — depend 
on the properties of the models. 

One aspect of major importance that has been considered 
by most works involving diagnostic tools to study the inner 
regions of stars is the ability of these tools to determine the 
stellar age. The correlations between the quantities studied in 
^ and the age of the models are clear. Hence, in principle 
the former might be used to infer the evolutionary status of 
pulsating star s. In fact, this co rrelation is responsible for the 
one found bv iMazumdar et aP (2006) between their quantity 
6 and the central hydrogen abundance Xc, (cf. figure 7 of their 
paper). Their 9 is proportional to the difference (D13-D02) 
averaged over several radial orders n, and through it reflects 
the change in the slope seen in the left panel of our Figure |4l 
Thus, our analysis also provide s the th eoretical basis for the 
correlation found by Mazumd ar et alj (j200Q). Nevertheless, 
it is important to investigate how these correlations depend 
on other stellar parameters, and to ascertain whether there are 
degeneracies in the age determination when the appropriate 
parameter space is explored. Also, as mentioned in § 13.11 
the models used to test our diagnostic tool did not include 
diffusion. The presence of diffusion will likely result in a 
smoother transition in c^. We hope to tackle the consequences 
of including diffusion in future work. 

Asteroseismic data with the precision necessary to test our 
diagnostic tool observationally are soon expected from the 
CoRoT satellite jBaglin et alj|2006b . The second scheduled 
5-month run will inclu de the star HP 49 933, which exhibits 
solar-like oscillations ("M osser et al J 1200 5) and has approxi- 
mately the required mass and age to reveal the expected sig- 
nal. Future observations of th e slightly more ma ssive spec- 
troscopic binary star 12 Bootis (iMiglio et al.ll2007h could also 
provide an interesting test of our predictions. 
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APPENDIX 

A. SOLUTION TO THE WAVE EQUATION NEAR THE INNER TURNING POINT 

In this appendix we derive an approximat e solution to E g. (|2]l, valid near the inner turning point z\ at which = 0. We closely 
follow the de rivation presented in §4.8.5 of lGoughl (1 19931) . though the dependent and independent variables in our Eq. (|2|i differ 
from those in lGoughl 's Eq. (4.8.5). 

We start by defining new independent and dependent variables, x and ?/; respectively 



sgn(K2) 



3gn(K2); 



\K\dz 



i^ = \x\-^\K\--S = s-2-~. 
With these definitions, x = at the inner turning point. Substituting these into Eq. ^ we find 



(Al) 
(A2) 

(A3) 
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where we keep the independent variable z in the right hand side (RHS) to keep the equation in a condensed form. 

To find an asymptotic solution to Eq. (IA3I) we regard the term on the RHS as a small perturbation to the Airy equation (lOlverl 
|1974|) . We note that our original transformati on of the inde pendent variable z = ln(r/ r,ef), was motivated precisely by the need to 
keep the RHS bounded when r tends to zero ( iLangeijri937h . In fact, a straightforward analysis of this term shows that it will tend 
to zero as r tends to zero. 

To leading order, we thus have that the solution to Eq. ( IA3I ) around the inner turning point can be obtained from 

^KiaAi (-x) + bBi i-x) , ( A4) 

where A, and B, are the solutions to the Airy equation, and a and b are constants. Since the solution must tend to zero as x tends 
to -oo, we must have b = 0. Using the ascending series for A, (e.g. Abram owitz & Stegun.l972t, §10.4, p.446), we thus find for 
\x\ < 1 

ipKa[Ai{Q)-A\(0)x+O{x^)], (A5) 

where 

and r is the Gamma function. Transforming back, we then find that 

S « a si [Ai (0) - A'i {Q)x] , (A6) 

with s defined by Eq. (IA2l l. 

B. VARIATIONAL ANALYSIS 

In this appendix we derive the expression for the frequency perturbations resulting from the presence of a small convective 
core. Our starting point is Eq. (Q, with 5I\ and Sh defined by Eqs. ([8]l and (|9]l. Introducing the new variables x and into 
Eqs. ^ and Q and noting that 

2 J j2^2 a2„i. a2. 



we find, after some algebra. 



dV- 
dx 



21] luSlo ! 



.^d^ and ^^-xV- 
2 dx^ dx^ dx^ ' 



dx dx^ 



dx, 



where Xq and Xr are the limits of the original integral in terms of the new variable x and 

21 



x + 



1 dq 
q dx 



Siip) [36{'fp)-4Sp]-u} r sSp . 

dx 



f 1 3 -1 
f2 = -r s — 
2 djc 



5{ip), 



h = -Ps-'qH(-fp), 



with q = s^l^r^l^r.^f'^l^p-^l^c-^ . 
Next, we express Eq. ( IB2l i in terms of the perturbations 5{^p)/"ip and 5c^ /c^. Using the relations 

d{jp) 5c^' 



Sp S{jp) Sc^ dSp 

— = J- and — = -pgrs 

p jp dx 



IP 



we find 



where 



5J'o = -{3r'q^-tp 
5Tx = 



Xr 
Xo 



dx 



diP^ 



dV 
1^ 



STi + — (5 Jo) H' + '5-^2 -7- + SJ^3 



dx 



dx, 



5{ip) 



IP 



x + 



1 dq 

q doc 



7/? - prs {\g + cj^r) 



+ 37/? 



da 



IP 



+ \r'(i'ps (4g + w^r) 



&2 



(Bl) 



(B2) 



(B3) 



(B4) 



6T2-- 
5Ty- 



1 3 -id?' 

-'^3-12 

-rs q -/p 



5{np) 

IP 

S(np) 
IP 
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The perturbation associated with the discontinuity at the edge of the convective core is rather locaHzed. Thus, we consider an 
interval [Xa,Xb] in which the perturbation differs from zero and assume it to be equal to zero outside that interval. By requiring 
both that the sound speed is everywhere the same in the perturbed and unperturbed models except near the discontinuity, and 
that both models are in hydrostatic equilibrium, the small difference in their density profiles extends beyond the region of the 
discontinuity, towards the center of the star Nevertheless, because we are interested in isolating the frequency perturbation 
associated with the rapid variation at the edge of the core, we do not consider the perturbation arising from this additional 
difference between the two models. 

Integrating by parts all terms of the integral on the RHS of Eq. ( IB4I ) except for that involving 6J-q, we find 



dx 



^2_|_d^d?// 
dx dx 



dx, 



(B5) 



where tp^ = J ip^dx. 

At the edge of the convective core, Sc^/c^ is discontinuous while all other perturbations in Eq. ( IB4I ) are not. Hence, the main 
contribution to the frequency perturbations will come from the term involving the perturbation to the sound speed. Accordingly, 
we neglect all other terms and find 



21] luSlo ' 



d 
dx 



n1 



i\P-dx. 



(B6) 



We model the sharp variation j at the edge of the convective core, xj, with a function of the type dc^ / oc [H(x—Xi)— 1 /2], 
where H is the Heaviside Step Function and the proportionality constant is determined by the amplitude of the sharp variation. 
Under this assumption, and keeping in mind that the term multiplying bc^ jc^ in Eq. ( IB6b is a slowly varying function of x, we 
have 

^A4rVpi(4g + c^2^)],5(x-Xd), (B7) 



djc 



rVps(4g + 



where A5 is the amplitude associated with the (5-function. We thus find, after some algebra. 



2/1 Lobuj ' 



-A5fl2A,.(0)2 



rrefc2|/:|2 



A,'(0) 2^ia;(0)2 

- — XX H ^ 

Am ' 3Ai(0)2 



(B8) 



In practice, if the approximations introduced above are reasonably satisfied, A5 is expected to be proportional to the amplitude of 
the sharp variation 5c^ /c^ at Xd, with the proportionality constant independent of the model considered. Indeed we found this to 
be the case in all of our models except for those older than ~4 Gyr. For these models the discontinuity in is not as sharp as for 
younger models. Even so, the functional form of the perturbations was still found to be relatively well represented by Eq. (IB8I 1 in 
this case. 

Finally, the constant a can be expressed in terms of the amplitude of the eigenfunction ^ at the inner turning point. Combining 
the definitions of t/j and S, we find 

m,=,,= [q-'(\^__^ = aAm- (B9) 
Substituting this expression for a into Eq. (IB8I 1. we finally have 



5ijJ : 



2ujI\ 



2(2 



1 1/2 



A-O 



A[{0) 2 lA[iOf 

w^-a;(0)W^^3a;(0)^ 



(BIO) 



For practical reasons, when calculating the theoretical frequency perturbations from Eq. ( IBIQI ). we use the quantities cj, Ii and 
i^A-o calculated for the model with the discontinuity. To first order, this should have no impact on the theoretical signal. 
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